Dipolar Ordering and Quantum Dynamics of Domain Walls in Mn-12 Acetate 
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We find that dipolar interactions favor ferromagnetic ordering of elongated crystals of Mni2 
Acetate below 0.8 K. Ordered crystals must possess domain walls. Motion of the wall corresponds 
to a moving front of Landau-Zener transitions between quantum spin levels. Structure and mobility 
of the wall are computed. The effect is robust with respect to inhomogeneous broadening and 
decoherence. 
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I. INTRODUCTION 

Molecular magnets exhibit quantum dynamics at the 
macroscopic level. The best-known expression of such a 
dynamics is the staircase magnetization curve that one 
observes on changing the magnetic fieldjii^i^ The steps 
occur due to Landau-Zener transitions between cross- 
ing quantum spin levels.** It has been previously demon- 
strated that dipole-dipole interactions in molecular mag- 
nets lead to ferro- or antiferromagnetic ordering of spins 
at low temperaturei^iii^i^ The Curie temperature as 
high as 0.9 K was reported in neutron scattering exper- 
iments on Mni2 Acetatei^ In this paper we re-examine 
the effect of dipolar interactions in Mni2 within numeri- 
cal model that treats spin-10 clusters as point magnetic 
dipoles located at the sites of a body centered tetrago- 
nal lattice. We find that elongated crystals must order 
ferromagnetically below 0.8 K. 

It has been noticed in the past^° that magnetic re- 
laxation in molecular magnets is a collective effect. In- 
deed, the change of the spin state of one molecule results 
in the change of the long-range dipolar field acting on 
other spins. When this change in the local dipolar field 
causes crossing of spin levels at a certain crystal site, the 
spin state of the molecule at that site changes as well. 
Quantum many-body Landau-Zener dynamics of molec- 
ular magnets has been intensively studied in recent years 
by means of Monte Carlo simulation o^^i*^'^'^ and by ana- 
lytical methodsii^ii^ In this paper we employ analytical 
model that takes into account both, local spin transitions 
and the long-range dynamics of the dipolar field. Within 
such a model it becomes obvious that existing Monte 
Carlo simulations of collective spin dynamics of molec- 
ular magnets have missed an essential feature of that 
process: Below ordering temperature the relaxation may 
occur via propagation of a domain wall (DW) separating 
spin-up and spin-down regions. Unlike domain-wall mo- 
tion in conventional ferromagnets, the dynamics of the 
domain wall in a molecular magnet is entirely quantum. 
It is driven by quantum transitions between spin levels 
that are crossed in a deterministic manner in space and 
time by a propagating wave of the dipolar magnetic field. 
Note that a propagating front of the magnetization re- 
versal has been recently observed in Mni2 crystals and 
interpreted as magnetic deflagratioufiSiiiii^ The latter is 



a classical phenomenon equivalent to the flame propa- 
gation, with the Zeeman energy playing the role of the 
chemical energy. Quantum mechanics enters the defla- 
gration problem only through the reduction of the energy 
barrier near the tunneling resonance. On the contrary, 
the phenomenon described in this paper has quantum ori- 
gin. It corresponds to a wave of Landau-Zener transitions 
generated by dipole-dipole interaction between magnetic 
molecules in a crystal. 

We consider quantum tunneling between two nearly 
degenerate ground states \iS) of magnetic molecules at 
low temperatures, interacting with each other as mag- 
netic dipoles. The molecules are then effectively de- 
scribed by spin 1/2 instead of spin 10. In the absence 
of quantum tunneling between |±5) these states do not 
communicate with each other so that any initial distri- 
bution of molecules in spin-up and spin-down sates will 
be preserved. The measure of communication between 
\±S) is their tunnel splitting A. The effects described 
in this paper, such as ferromagnetic ordering and motion 
of domain walls, can be observed only if A is sufficiently 
large. In Mni2 it can be controlled by the transverse 
magnetic field. Since we are interested in the motion of 
domain walls, we choose elongated sample in the shape 
of a long cylinder of length L and radius i?, the quan- 
tization axis of spins being directed along the z-axis of 
the cylinder. We restrict our consideration to the states 
only weakly nonuniform at the lattice scale, so that spins 
in macroscopic regions are parallel to each other. This 
can be achieved by either polarizing spins by the exter- 
nal magnetic field or through ferromagnetic order which, 
as we shall see below, plays an important role al low 
temperatures where many experiments were performed. 
We further simplify the problem by ignoring inhomo- 
geneities along the perpendicular axes x and y, so that 
az = (Sz) /S depends on z only. 



II. THE MODEL 
A. The density matrix equation 

The effective Hamiltonian of one magnetic molecule 
at site i within the MFA, taking into account only the 
two ground states \±S) of a molecular magnet at low 
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temperature, can be formulated in terms of pseudospin 
(below spin) a as 



The equilibrium solution of Eq. ([S]) has the form 



(1) 



Here W — 2Sgfi^B2 is the energy bias that generally de- 
pends on time via the total longitudinal field including 
the external and dipolar fields. A is tunnel splitting de- 
fined by the uniaxial anisotropy D and the terms in the 
Hamiltonian that cause tunneling, e.g., transverse field 
B±. Finally, ct^, ax are Pauli matrices. The energy levels 
of this Hamiltonian for an instantaneous value of W are 



A2, 



(2) 



where loq is the corresponding transition frequency. 

The density-matrix equation (DME) for the spin in the 
time-dependent adiabatic basis, formed by the instanta- 
neous eigenstates \x±) of HcS, has the form 

^P++ = ((x+ |x+) + (x+ |x+)) P++ + 
(x+ \x-)p-+ + p+-{x- |x+) - ^-+p++ + r+-p__ 
^lP+^ = ((x+ |x+) + (x- |x-)) P+- + p++{x+ |x-) 



dr 

+(x+ \x-)p-- 



(3) 



where F |_,F_| are up and down relaxation rates for 

the levels e±, satisfying the detailed balance condition 
F+_ = g-?ic^o/(fcBT)p_^ rpj^g elements of the density ma- 
trix satisfy + p__ = 1 and p ^ = (p^ ) . 

Taking time derivative of the spin expectation value 
<T = Tt(p&) and using the relation 

cr = {p^_ + p_^) e^+i {p+- - P-+) ey+{p__ - e^, 

(4) 

one finds that in the chosen time-dependent frame the 
DME describes damped precession of cr about the effec- 
tive field u3q + Osy, where 



h 



(5) 



and cos 6* — W/ \JW^ -I- A^ describes the orientation of 
ujq. Switching to the laboratory coordinate frame (which 
amounts to dropping non-adiabatic term d in the effec- 
tive field) one obtains 



CTq = o-Q- 



fjQ = tanh 



(7) 



Remember that the equations above are written for the 
spin on a site i, although the index i is not explicitly 
written for brevity. The equation of motion ([6]) for spin i 
couples to those for all other spins j via the magnetostatic 
contribution in W considered in the next section. 

In transverse field B\^ satisfying hwQ <ti gpQB± the 
relaxation rate F due to direct processes is given by 



F = 



UirEf 

>l/4 



2fcBT' 



(8) 



where = {^pv^h^^ is a characteristic energy. For 
Mni2 Et/kB ^ 150 K but the factor 1/2 should be in- 
troduced in F since there is only one transverse sound 
mode. Usually direct rates are proportional to cjq but in 
the case of tunneling there is an additional factor ljq ^ 
accounting for the decrease of hybridization of the two 
levels with increasing the energy bias W. 

The model formulated above describes the Landau- 
Zener effect in the case of time-dependent energy bias 
W that crosses zero, as well as relaxation. In the pres- 
ence of dipolar coupling this model describes magnetic 
ordering and domain-wall dynamics. Especially inter- 
esting is possible interplay between the DW dynamics 
and LZ effect that can possibly lead to Landau-Zener 
fronts. The model can be extended by including heating 
the sample as the result of spin relaxation. This leads 
to formation of deflagration fronts if the relaxation rate 
strongly increases with temperature. 



B. The dipolar field 

The energy bias W in the equations above at the site 
i is given by 

= 2SgpB {Bz + Bf^^ ) = VFe,t + W^""^ , (9) 

where B^ is the z component of the external field and 
Bf'J is the dipolar field at site i. The dipolar component 
of the bias is given by 



(T = [cr X Wg] 

F / Wo ■ <T 

— — I (T — 

2 



r — —^a 



(To 



(6) 



where F = F |_-|- F^ and ctq is the thermal equilibrium 

value of the pseudospin, corresponding to the instanta- 
neous value and direction of ujq. The second term in this 
equation corresponds to the relaxation of the spin com- 
ponent perpendicular to ljq while the third term corre- 
sponds to relaxation along ljq, the latter being twice as 
fast as the former. 



where E^ = (gp^S) /vq is the dipolar energy, vq is the 
unit-cell volume, and 

is the dimensionless dipole-dipole interaction between the 
spins at sites i and j ^ i. The z component of the dipolar 
field itself is given by 



^(D) ^ SgpB j^ 
vo ' 



(12) 
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To calculate the dipolar field, one can introduce a 

1 /3 

macroscopic sphere of radius ro satisfying Vq rg <^ L 
aroimd the site i, where L is the (macrocopic) linear size 
of the sample. The field from the spins at sites j inside 
this sphere can be calculated by direct summation over 
the lattice, whereas the field from the spins outside the 
sphere can be obtained by integration. The details are 
given in the Appendix. In particular, for a uniformly 
magnetized ellipsoid the total result has the form 



Dzz 
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independently of i, where 



Dzz = ^iT''^ + 4vrj. (l/3 



(13) 



(14) 



and v is the number of molecules per unit cell. For the 
demagnetizing factor one has n(^) = 0, 1/3, and 1 for 
a cylinder, sphere, and disc, respectively. One obtains 
D^zz^^ = for a simple cubic lattice, < for a 

tetragonal lattice with a = h > c, and D'tF^^ > for 



b < c. 



Eo = ^{1/2)DzzEd 



(15) 



is the dipolar energy per site. 

The two best known molecular magnets are Mni2 and 
Fes, both having total spin S — 10. Mni2 crystal- 
lizes in a body-centered tetragonal (bet) lattice with 
a = b = 17.319 A, c = 12.388 A (c being the easy axis) 
and the unit-cell volume — obc — 3716 A^, with two 
molecules per unit cell, v — 2. Fes has a triclinic lattice 
with a — 10.52 A (a being the easy axis), b = 14.05 A, 
c = 15.00 A, a = 89.9°, /3 = 109.6°, 7 = 109.3° and 
Wo = a6c sin a sin /3 sin 7 — 1971 A'^. The characteristic 
dipolar energies thus are E-o/k-Q — 0.0671 K for Mni2 
and E^/k^ = 0.126 K for Fes. 

For Fes direct numerical calculation yielda^^ o'tf^'^ — 
4.072, thus for the cyhnder Eq. HH) yields Z)i7'^ = 
8.261. Our result Eq = -^^.UIEd for the elongated Fes 
crystal is in qualitative accord with Eq = —A.IQEd of 
Ref. i 

For Mni2 one obtains D^zz^^ = 2.155 that results in 
= 10.53 for a cyhnder. Then Eq. HID yields the 

dipolar field ^ 0.0526 T in an elongated sample. 

On the top of it, there is a weak ferromagnetic exchange 
interaction between the neighboring Mni2 molecules that 
creates an effective field 7 G from each neighbor With 
8 nearest neighbors in the bet lattice, this effectively adds 
1.12 to Dzz in the ferromagnetic state. Thus for Mni2 
one obtains effectively o'fz^^ ~ 11.65. 

It can be shown that all other types of ordering have a 
lower value of Dzz and thus a higher value of the ground- 
state energy Eq for both Mni2 and Fes of a cylindric 
shape. For Fes the state with ferromagnetically ordered 
planes alternating in the c direction has Dzz = 8.18, the 
dipolar field being shape independent. This value is very 
close to 8.261 for the ferromagnetically ordered cylinder 



and it has only slightly higher energy Eq. The states with 
ferromagnetic planes alternating in the b direction has 
Dzz — 8.12, while the state of ferromagnetic chains di- 
rected along the a axis and alternating in b and c direction 
has Dzz = 8.01. One can see that a crystal of Fes cooled 
in zero field will show a random mixture of different types 
of ordering. On the other hand, as different kinds of or- 
dered states are separated by energy barriers, prepared 
ordered states should be robust metastable states. In 
particular, ferromagnetically ordered state of an elon- 
gated Fes crystal, obtained by cooling in the magtic field, 
should be stable after removing this field at low temper- 
atures. 

For Mni2, states with ferromagnetically ordered planes 
alternating in the a or 6 directions in each sublattice have 
Dzz = 9.480, independently of the shape and of the ex- 
change interaction. The state with alternating chains 
directed along the c direction has a very close value 
Dzz — 9.475. For the two-sublattice antiferromagnetic 
ordering one obtains Dzz — 8.102. All these values are 
essentially lower than the dipolar field Di'^/^^ ~ 11.65 for 
a cylinder. Thus one can expect that elongated crystals 
of Mni2 will order ferromagnetically without competition 
of other states. 

On the other hand, for a spherical shape the states with 
alternating ferromagnetically ordered planes and chains 
are energetically more favorable for both Fes and Mni2. 

For a cylinder magnetized with cr^ = crz(z), the field 
along the symmetry axis has the form 

Dzziz) = ^ dz'—^^^^^^^^-kaz{z), (16) 



[z'-zf+R^ 



where 



k = 87ri//3 - D'iF'^^ = Attv - Li^^'^ > 0, (17) 



k — 14.6 for Mni2 and k — 4.31 for Fes. For a uniformly 
polarized cylinder Eq. (fT6)) yields 



Dzz{z) = 2t:v 



z + L/2 



^^{z + L/2f 



i?2 



L/2 



(z-L/2f + R\ 



Gz - kaz{\%) 



In the depth of a long cylinder, L':^ R, the dipolar field 
is Dzz = (^oi'^z^^ +AT:i^/3j az, in accordance with Eq. 

At an end of a long cylinder one has Dzz — 
(^D'fF^^^ — 2ttv/3^ Uz that can have the sign opposite to 
that in the depth of the sample. In particular, for Fes 
with b^zf''' = 4.072 and 1/ = 1 one has Dzz = IMuz- 
Thus a homogeneously magnetized state in zero field is 
stable. To the contrary, for Mni2 with o'fz^^ = 2.155 
and V — 2 one has Dzz = —2.Q3az- This means that a 
homogeneously magnetized state in zero external field is 
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unstable with respect to domain formation, beginning 
in the vicinity of the ends of the cyhnder. At some 
point near the end of the crystal the resonance condi- 
tion Dzz{z) = is satisfied that leads to spin tmmcling 
and the decay of the initially homogeneously magnetized 
state. 

If there is a domain wall at z = (cz — > tioo for z 
—oo) in a long cylinder that is narrow in comparizon to 
R, Eq. (fT6|) far from the ends yields 



-Ajtv 



zo 



^{z-z^f + R^ 



ka,{z). (19) 



Since the coefficient of the local term is negative, the lat- 
ter changes in the opposite direction with respect to the 
first term. This creates three zeros of D^ziz). However, 
such exotic DWs do not exist, as we will see below. Ther- 
modynamically stable domain walls have only one zero of 
Dzz{z)., thus their width at low temperatures is of order 
R. 



III. DIPOLAR ORDERING AND STATIC 
DOMAIN WALL PROFILE 

A. Ferromagnetic ordering 

Dipolar field causes ferromagnetic ordering in elon- 
gated crystals of molecular magnets that is described 
within the MFA by the Curie- Weiss equation following 
from Eqs. (HI), ([9]), and p3)) . This equation for a homo- 
geneously magnetized sample of an ellipsoidal shape has 
the form 



We. 



2Ey)Dzz(Tz ^ , ^0 
— tanh ■ 



2kBT' 



(20) 



where Hujq = \J {W^^^t + '^Ej^Dzztr z) +A2. This equa- 
tion is similar to that for the Ising model in a transverse 
field, here A. The Curie temperature satisfies the equa- 
tion 



1 



2Et,Dz 



■ tanh ■ 



2kBTc 



(21) 



This equation has a solution only for A < E^Dzz as the 
transverse field tends to suppress the phase transition. 
In the actual case A <C E-qDzz one obtains 



Tc = EoDzz/kB- 



(22) 



For a Mni2 one has Ej:, /fce — 0.0671 K. Thus for a cyhn- 
der {Dzz ^ 11.65) Eq. (HH) yields Tq ~ 0.782 K. This 
is close to the value 0.9 K reported in Ref. Hi Note 
that the MFA does not account for fluctuations that usu- 
ally lower Tc- However, as we have seen above, for the 
Mni2 cylinder the main contribution to the dipolar field 
comes from a great number of distant spins, thus the 
MFA should work well. The small disagreement in Tc 
can be attributed to approximating magnetic molecules 
by point dipoles in our calculations. 



The equilibrium value of cr^ in zero field below Tq sat- 
isfies the equation 



2EuDzz ^ , fuvo 
1 = — : tanh ■ 



hun 



2kBT' 



(23) 



where hujo = y {2EjyDzzO'z) + A^. In the realistic case 

A ^ E-qDzz not too close to the Curie point one can 
neglect A^ in fiiOQ, then Eq. defining the spin polar- 
ization below Tc takes the form 



Gz — tanh 



kjiT 



(24) 



For Mni2 the ratio A/ (2EuDzz) is small. Even for B± = 
5 T that corresponds to h± = g^^B±/{2SD) — 0.61, 
one has A/kn =0.00166 K and thus A/ {2EbDzz) ^ 
2.5 X IQ-^. 



B. Static domain wall profile 

The magnetization profile (Tz{z) in a domain wall join- 
ing regions with dz = ± (Too = ±0'^'^ in zero field satisfies 
the equation 



2EbDzz{z) ^ , hujQ 
a, — tanh ■ 



2fcBT' 



(25) 



where Hloq = y [2EB)Dzziz)] + A^. Here Dzz{z) is given 
by Eq. (fT6|) that makes Eq. ([25|l an integral equation. 
By symmetry, Dzz{z) goes through zero in the center 
of the DW. At very low temperatures, T <C A/fce, the 
argument of tanh is always large, thus Eq. (|25p simplifies 
to 



2E^Dzz{z) 



[2E^Dzz{z)Y + 



(26) 



In this limit the spin length is nearly constant, cr = 1, so 
that Ox — — c^. The spin-polarization deficit in the 
wall is defined by 



1 — cr = 1 — tanh 



{2E^Dzz{z)X + A2 



2fcBT 



2fcBT 



{2E^Dzz{z)X + 



< 1 (27) 



and it reaches the maximum at the DW center. With the 
help of Eq. one finds 



\-u = ——^\-al. 



(28) 



In the opposite limit Ajk^, ^ T Eq. (pS)) simplifies to 

ExiDzzXz) 



Oz — tanh ■ 



A:bT 



(29) 
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FIG. 1: Magnetization profile of a domain-wall in a Mni2 
cylinder at two diflterent temperatures. 
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Domain wall profile 
in a long Mn^2 cylinder: 
Dipolar field 
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FIG. 2: Dipolar-field profile of a domain-wall in a Mni2 cylin- 
der at two different temperatures. 



The transverse spin component in this case is given by 



tanh 



(30) 



and cr^ = A/(2A:bT) < 1 at the center of the DW. For 
Mni2 is very small except of very low temperatures 
and very large transverse fields. Such a domain wall is 
the linear (Ising-like) domain wall found in ferromagnets 
in a narrow temperature range below Tq ^^0i^1i22i23 ^ ^qYI 
as in the low-temperature strong-anisotropy ferromagnet 
GdClsj^ Absence of a strong short-range exchange inter- 
action in molecular magnets makes domain walls linear 
practically in the whole range below Tq. 

For Mni2 the splitting A is typically so small that one 
can consider the limit T ^ of Eq. (p9| . Then three 
regions arise: (i) left from the DW where az — 1 and 
the argument of tanh is infinite, that is, Dzz{z) > 0; (ii) 
inside the DW where \az\ < 1 and thus Dzz{z) = 0; right 
from the DW where Cz = — 1 and Dzz{z) < 0. Using Eq. 
(fTB)) with L = oo one obtains the integral equation for 



the nontrivial region inside the DW centered at z = 0: 
27rR^crz{z') 



= 1/ 



dz'- 



{z' - zf + i?2 



I 



3/2 



k(Jz{z) 



z + / 



^{z-lf^R^ ^{z^lf^B? 



,(31) 



where / is the DW width that should follow from 
the equation above and the boundary conditions are 
(Tz(±Z) — ^1. One can see that still in for T ^ the 
DW profile results from an essentially integral equation 
that does not have a general analytical solution, although 
analytical solutions in limiting cases do exist. Rewriting 
Eq. (pTj) in terms of i = zjR^ one can show that the DW 
profiles scales with the cylinder radius i?, i.e., / ^ R. 

Numerical solution using relaxation to the equilibrium 
described by Eq. (j25p with proper boundary conditions 
at the ends of a long Mni2 cylinder is shown in Figs. [1] 
and [21 One can see that at T = the magnetization 
profile turns to ±1 beyond the region of the DW of the 
length / ~ i?, in accordance with Eq. ([51]) . At finite 
temperatures the solution for (Tz{z) is more resembling a 
tanh. The DW width I increases with temperature and 
diverges at Tq- The solution for the dipolar field in Fig. 
[2] shows that Dzz{z) ^ at T ^ inside the domain 
wall while it is approaching the asymptotic value ±11.65 
in the domains. At finite temperatures the dependence 
Dzz{z) smoothens out. Overall the region of inhomo- 
heneity of the dipolar field is broader than that of the 
magnetization. 

The DW width I can be defined as a slope 



1 daz 

Too dz 



(32) 



z=0 



Some results for / can be obtained analytically. At T <C 
Tc , assuming that the magnetization profile in the DW is 
close to the piece-wise linear, one obtains from Eq. 

r X -2 n ^1/2 

n2 / " 

-R ^ [-E^ + \ 

_ AniyEu + T - Tc 



(33) 



in a qualitative agreement with numerical results in Fig. 
H At T Tc the solution of Eq. ^ is 



(Jz{z) = (Too tanh 



where CToo — V3{Tc/T — 1)^/^ and /ht satisfies 



Z?rp In ■ 



HT 



with /lt of Eq. 



diverging as 



^LT ^ / 27r^D 

R "V^B(rc-T)- 



(34) 



(35) 



(36) 
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Domain wall width / 
in a Mni2 cylinder of raduis R 




FIG. 3: Temperature dependence of the DW width / and its 
low- and high-temperature forms in a Mni2 cylinder. 



Temperature dependence of the approximations Zlt and 
Zht, together with numerical result for the domain wall 
width /, is shown in Fig. [31 



IV. DOMAIN- WALL DYNAMICS 

In the sequel we will be mostly interested in the motion 
of domain walls with a small speed, induced by a small 
external bias Bz- In this case the system does not deviate 
much from the equilibrium. In this case domain walls in 
a sufficiently long sample are stationarily moving with a 
speed Vuw that is proportional to the bias field. This 
relation has the form 



(37) 



where /i^w is the linear mobility of the DW. At high 
values of the dependence vuwiHz) becomes nonlin- 
ear. The linear mobility /iD-^ can be calculated using the 
static DW profile discussed above with the help of the en- 
ergy balance argument. To this end, the LLB equation 
should be transformed into a special form near the 
equilibrium. 



A. Landau-Lifshitz-Bloch equation near 
equilibrium 

One can project relaxation terms in Eq. © on the 
directions parallel and perpendicular to cr using 



R 



(cr-R)cr [crx [cr x R]] 



(38) 



cr-' 

This results in a Landau-Lifshitz-Bloch (LLB) equation 
& — [cr X LJo] 

+r 



i^o ■ cr /CTo _i\ (i_ '^o • 
ujo<j \ (T / 2 \ ujo<j 

ljq ■ cr\ [crx [cr X ljo]] 



(39) 



Here, in contrast to ferromagnets where cr = 1 is enforced 
by a strong exchange, a can essentially deviate from 1 
because of the first relaxation term. Close to equilibrium 
vectors ujq and cr are nearly coUinear, thus • cr = LUoa 
holds up to quadratic terms in small deviations, whereas 
cr is close to (Jo- Thus the equation above simplifies to 

/o-Q \ F [crx [cr X ujn]] , , 

d-= [crxa;o]+F(^-l cr--i ^ ^. (40) 

\ cr / 2 aujo 

Note that ctq is formally defined by Eq. ([7]) and it is a 
function of uq that can depend on time. Thus, in gen- 
eral, CTq corresponds not to a true equilibrium but to 
an instantaneous equilibrium to which the system tends. 
If the state of the system, as well as the instantaneous 
equilibrium, is close to a true equilibrium, one can ap- 
proximately express ujq through cr and small & that has 
an important application. First, multiplying Eq. (1401) by 
cr one obtains the relation 



CTo = CT + 



l {<T-cr) 
r cr ' 



(41) 



Second, using ivq ■ cr = luqct one can rewrite Eq. (|40p in 
the form 

(T = [cr X Wo] +r (—- l) cr-^Jcr+ J— a;o. (42) 
\ cr / 2 2 UJQ 



From the latter follows 

[cr X <t] = crujQa — ujQa 



2'^-[^xa;o]. (43) 
2 ujQ 



Eliminating [cr x ujq] and using Eq. (j4ip . one obtains 



cr X cr 



i^a = t^o — 



^or/2 



crx cr X cr 



(44) 



Here the scalar ojo can be found from Eqs. ([7]) and ((4T 
2kBT 



h 

2fcBT 

h 

2kBT 



arctanh (erg) 

/ 1 (cr • d") 

arctanh cr H 

V F cr 

2fcBT (T 

arctanh (cr) H -■ (45) 

nV 1 — cr^ 



The & correction here should be taken into account in 
the first term of Eq. (|44)) . 



B. Linear mobility of domain walls 

The time derivative of the magnetic energy of the sam- 
ple per unit cross-sectional area due to dissipation has the 
form 



U 



dz hijJo{z) ■ cr{z). 



(46) 



On the other hand, motion of the DW results in the 
change of the energy of the spins in the external field 
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at the rate —a^WexiVDW^ where (Too > is the spin 
polarization in domains and Wext is given by Eq. ([9]). 
Equating this rate to Eq. (pS]) . one obtains the energy 
balance relation 



CooWcxt^DW = 



dz hijJo{z) ■ <t{z). 



(47) 



This relation allows one to obtain the linear mobility of 
a domain wall as an integral over the static DW profile 
without solving a complicated problem of dynamical cor- 
rections to this profile. To find the linear mobility, one 
has to express uJo{z) through <t(z) using Eq. (jU]) that 
yields 



cr ■ cr cr^o"^— (cr • or) 



: UJQ- 



r2/4 



(48) 



After that one has to use the fact that for a domain 
wall stationarily moving in the positive z direction, all 
quantities depend on the combined space-like argument 
^ = z — vuwt, so that 



-vuwdcr/dz. 



(49) 



Using Eq. ([i5|) . one obtains the energy balance equation 
in the form 



fcxjWoxt 

dz I ^2A:Brarctanh ((t) 

arctanh (cr) 
1 / do- 



er 

a dz 

1 
a 



r 1 - (t2 rfl 

AV4 
arctanh^ (cr) 



da 

dz 



where 



A 



dz 



(50) 



(51) 



Here the first term in the rhs can be easily integrated and 
gives a zero contribution. Thus the rhs of Eq. (|T7)) be- 
comes proportional to v^^. This is very fortunate since 
now one does not have to take into account deviations 
from the static DW profile. In all practical cases the in- 
equality A ^ 1 is strongly satisfied, thus one can replace 
[. . .] ^ ^ 1. After that for the DW speed in the linear 
regime one obtains 



Sa 

OO 



2fcBT k^T 
where is the characteristic speed defined by 



(52) 



dz- 



1 



1 



r I 1 - 
AV2 



cr arctanh (cr) 



da 

dz 

da 

dz 



da 

dz 



,(53) 



that resembles the expression for 

ferromagnetsj22i2ii2ii2^i2S Note that T is not a con- 
stant and it has to be kept in the integrand. One can see 
that if cr —> 1 in the domain wall, as is usually supposed 
to be the case in ferromagnets, then da/dz ~^ 0, only 
the second term makes a contribution, and /^qw oc l/F. 
However, if a only slightly deviates from 1, then because 
of the very small A the second term becomes irrelevant 
and one obtains /i^w oc F from the first term. For the 
linear (Ising-like) domain wall, a ~ aze^, the expression 
in the square brackets in second term of Eq. (j53p 
disappears. This is the case for Mni2 except for very 
large transverse fields and very low temperatures. 

The crossover between the two different types of the 
mobility behavior occurs in the low-temperature range 
T <C A/Zcb, where the deviation of a from 1 is small and 
given by Eq. ((28)) . Here in Eq. ([53]) one can simplify 



1-a'' 



da 
dl 

da 

dz 

arctanh (cr) 



2(1 



4fcBr 

A 

da^ 

dz 



4fcBr 

A 

„2 



'da^ 
dz 



do_z 

dz 



dz 



im- 

2 1 



- a 
A 



ii„(_ 

2 \2kB.T 



-1/2 



(54) 



This yields 



dz 



1 4fcBT al 



A H-alf" 



' da^ 

dz 



(55) 



One can see that the second term dominates and the DW 
mobility has the form /Iq^y o<^ 1/F only at extremely low 
temperatures. 



T < 



2/3 



(56) 



up to the log term. In most region one has /Idvv °^ ^■ 

For Mni2 in the typical case A <C the DW is Ising- 
like, and one should use the first term of Eq. ([53]) . where 
ax and ay are negligibly small, so that in Eq. (|5^ one 
has 



dz- 



1 1 



Fl-cr? 



da^ 
dz 



(57) 



In the case of the unspecified relaxation rate F that is 
considered as a constant instead of Eq. (|52p one obtains 



kjiT 



(58) 
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with 



dz 



1 [da. 



(59) 



At T <C Tc according to Eq. (|3T|l one has /* ~ i? inde- 
pendent of Ed with a coefficient dependent only on the 
lattice structure. Thus obtains /^dw ^ 1/^ low tem- 
peratures. In this case numerical calculation yields 
diverging both at T — > and T — > Tc and having a min- 
imum at about 0.25Tc. One can improve the numerical 
procedure eliminating division 0/0 in domains at T ^ 
using Eq. ^ and 



— — arctanh (cr^) . 



1 — a1 dz dz 



This yields 



Ed 



dz 



doz dDzz 
dz dz 



(60) 



(61) 



This expression is very convenient since at low tempera- 
tures Dzz ^ k-oT/ED inside the DW that makes no com- 
putational problems. The non-monotonic dependence of 
wdw on temperature should be taken with caution, how- 
ever, because it is related to our assumption that F is a 
constant. 

For F given by Eq. ([8]) is intimately interwoven with 
the rest of the mobility formula. Using Eq. in the 
form 



^(^).tanh.^-" 



hwo 2kBT' 
one can transform Eq. ([5]) into the form 



F = 



127rEf 



(62) 



(63) 



hcTziz) ' 

eliminating W{z) with the use of Eq. one obtains 



S^A^ [g^BB^Y SfceT arctanh [uz] 



UnEf 



Now Eq. (|52p can be rewritten as 



where 



r, -p(cxt) , 



12nEf 



(64) 



(65) 



(66) 



is the prefactor in Eq. ([5]) with fkuo => Wcxt and is a 
characteristic length defined by 



dz- 



da-, 



arctanh (cr^) 1 — ct^ \ dz 



(67) 



Similarly to Eq. ((57|) with F ^ 1, this integral is finite 
at T — !■ 0. To improve its numerical convergence at low 



I'Yy^lB, (relative units) 




0.8 j/r^LO 



FIG. 4; Temperature dependences of the DW mobility in the 
cases of F = const and of F due to the direct phonon processes. 



temperatures, one can rewrite it as 



dz 



a z da z darctanh [a Y) 
arctanh ((T^) dz dz 



<Tz dcTz dPzz .gg. 

Dzz dz dz 



Again, at low temperatures ^ R with a coefficient 
depending only on the lattice structure, c.f. Eq. (pT|) . 

Numerical results for temperature dependence of the 
velocity for F = const and for F given by a direct phonon 
process in a transverse magnetic field^i are shown in Fig. 

SI 

More generally, one can consider the relaxation rate of 
the form 



F = A{fiajoT coth 



hujQ 

2fcBT' 



(69) 



where A — const at low temperatures. In the case A <C 
Ed one has hujo ~ W{z), then with the help of Eq. 
one can rewrite F as 



F = A (2/cbT arctanh (cr^))" — . 



DISCUSSION 



(70) 



We have demonstrated that elongated crystal of mag- 
netic dipoles arranged in a body centered tetragonal lat- 
tice should exhibit ferromagnetic ordering. For Mni2 
Acetate the corresponding Curie temperature computed 
within mean-field model is about 0.8 K, which is close to 
the ordering temperature 0.9 K reported in recent neu- 
tron scattering experimentsi^ Other magnetic phases are 
separated from the ferromagnetic phase by a large energy 
gap, thus making ferromagnetism in Mni2 not difficult to 
observe. This suggests that the analysis of all previous 
data on Mni2 Acetate that were taken below 0.8 K should 
be re-examined for possible effects of ferromagnetic order. 

Dipolar-ordered crystals of molecular magnets should 
possess domain walls. In a long crystal of length L and 
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radius R <^ L, the typical width of the domain wall 
is of order R. When a small bias magnetic field 
is applied, the domain wall moves at a speed wdw ^ 
[Sg^i-gBz/ {k'QT)\(r)R, where (F) is the average spin re- 
laxation rate. At, e.g., 5" = 10, = 0.1 T, T = 1 K, 
and i? = 1 mm, this gives udw ^ 1 ni/s for (F) = 10'^ 
s-i and WDw ~ 10^ m/s for (F) = 10^ s'^. It should 
be emphasized that contrary to superradiance and laser 
effects in molecular magnets^i^ quantum dynamics of 
domain walls is robust with respect to inhomogeneous 
broadening of spin levels and phase decoherence of spin 
states. Crossing of spin levels due to a moving front of 
dipolar field is sufficient for the effect to exist. Neither 
very narrow spin levels nor phase coherence of spins in 
the domain wall are required. We, therefore, believe that 
this effect should not be difficult to observe in molecular 
magnets. 

Our conclusion that low-temperature ferromagnetic 
phase is unstable against division into domains may seem 
to contradict experiments on Mni2 in which, at low tem- 
perature, the crystal was shown to maintain finite mag- 
netization in a zero field for a very long time. This ap- 
parent contradiction becomes resolved if one notices that 
in the zero-temperature limit any evolution of the mag- 
netic state of the Mni2 crystal can only occur through 
quantum tunneling between spin levels. The measure of 
this tunneling is the splitting A which for low lying spin- 
levels is negligibly small in zero field. To increase A one 
should apply transverse magnetic field. Note that some 
evidence that transverse field is needed to achieve local 
ferromagnetic order was obtained in Ref. 9. In the ab- 
sence of the longitudinal field, the existence of the global 
ferromagnetic order would be hindered by the presence 
of domain walls. One would need local measurements of 
the sample to directly confirm ferromagnetic order below 
the Curie temperature. However, in the presence of suffi- 
cient transverse field, a relatively weak longitudinal field 
of the order of the dipolar field should be able to drive the 
domain walls out, thus resulting in the uniform magne- 
tization. Another simple experiment can be conducted 
to observe the domain wall entering the crystal. First 
the crystal should be uniformly magnetized by a strong 
longitudinal field. Than the longitudinal field should be 
switched off and a sufficient transverse field should be 
turned on. Finally, the compensating longitudinal field 
should be applied. As soon as the resonant condition is 
achieved, the domain wall should enter the crystal, re- 
sulting in the drastic change of the magnetization. Some 
of these effects must have been already observed in Mni2 
crystals but attributed to resonant quantum flipping of 
magnetic molecules at random sites, instead of coherent 
motion of domain walls. 

The motion of a domain wall should not be confused 
with magnetic deflagration . ^^i^^i^^ The latter is also char- 
acterized by a moving front that separates regions with 
opposite magnetization. However, while deflagration is 
driven by the release of heat due to a strong field bias, 
the motion of a domain wall is driven by Landau-Zener 
transitions in a relatively small bias. In experiment, the 
two effects can be easily distinguished from each other if 



one notices that the mobility of the domain wall depends 
strongly on the transverse field while deflagration must 
have little dependence on the transverse field. 

In the Introduction it was pointed out that domain 
walls in elongated crystals of molecular magnets are prop- 
agating waves of Landau-Zener transitions. The tunnel- 
ing and relaxational dynamics of this process is described 
by the density-matrix equation ([3]) that has a form of the 
equation of motion for the magnetization, Eq. ^ in the 
two-state lis*) approximation. The energy sweep in the 
LZ effect is due to the change of the dipolar field cre- 
ated by quantum spin transitions \S) ^ \~S) in the 
moving front. Thus, in contrast to the standard LZ 
effect with externally controlled time-linear sweep, the 
LZ effect considred here includes a self- consistent sweep 
that turns out to be time nonlinear. Indeed, the spa- 
tial profile of the dipolar field in the DW in Fig. [5] is 
strongly nonlinear at low temperatures. As the domain 
wall moves, a time-nonlinear sweep is generated. LZ ef- 
fect with time-nonlinear sweeps generated by interaction 
between two-level systems have been studied in Ref. HO- 
It was shown that if the time dependence W{t) becomes 
steeper near the resonance, the staying probability P in- 
creases, whereas in the opposite case the transition prob- 
ability 1 — P increases. Sweeps that become flat near the 
resonance lead to a nearly complete transition, see, e.g.. 
Fig. 7 of Ref. i30|. This is also the case for the sweep 
corresponding to the T ^ curve in Fig. [21 In this 
case the LZ transition is complete, P = 0, that is exactly 
what should be for a moving DW. Externally controlled 
sweeps leading to a complete transition have been studied 
in Ref. [U. The case of moving domain walls is another 
realization of the full-transition Landau-Zener effect. 

Finally we would like to comment on the square root 
of time magnetic relaxation observed in Mn-12 at low 
temperature. Existing theoretical works attempted to 
explain this behavior by flips of molecular spins at ran- 
dom sites of the crystal. iiii^ii^ They missed the fact that 
the low-temperature phase of Mni2 Acetate is ferromag- 
netic. Relaxation in the ferromagnetic phase must be 
determined by the motion of domain walls. Within such 
picture the square root relaxation may have simple ex- 
planation through random walk of the domain wall. If 
the wall moves randomly at a constant rate to the right 
or to the left, the resulting displacement and, thus, the 
change in the magnetization should be proportional to 
the square root of the number of steps, that is, to the 
square root of time. Such random walk may be thermal 
or it may be caused by the adjustment of the domain 
wall to local solvent disorder or to local hyperfine fields. 
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APPENDIX A: CALCULATION OF THE 
DIPOLAR FIELD IN A CRYSTAL LATTICE 

The dipolar field bI^-* created on the lattice site i by 
all other spins, j ^ i, aligned along the z axis can be 
written in the form 



with the dimensionless vector Di z given by 



3 



(Al) 



(A2) 



Here cr is the Pauli matrix and 



0, 



vo- 



n., ^ ^ (A3) 



with Vo being the unit-cell volume. As said above, Di z 
of Eq. (jAl|) can be represented as a sum of the contri- 
butions from the molecules inside and outside the sphere 

1/3 

ro around the site i satisfying <C ^ L, where L is 
the (macrocopic) linear size of the sample. The field from 
the spins at sites j inside this sphere rp can be calculated 
by direct summation over the lattice, whereas the field 
from the spins outside the sphere can be obtained by in- 
tegration. Replacing the index i by the argument r, one 
can write 



where 



and 



D^(r) =Di"P'^)(r)+D^(r) 



(A4) 
(A5) 



DUr) = -/ dV^(r-r')a,(r'). (A6) 

"^0 J|r-r'|>ro 

In Eq. (jA5|) it is assumed that the change of ct^ (r) inside 
the sphere tq is negligible. D^z^^^ is a constant depending 
on the lattice structure. In Eq. (jA6[) v is the number of 
molecules per unit cell. 



(r) = Vo 



3nr(e3 • rir) - 



= Wo rot 



[e^ X rir] 



(A7) 



and rir = r/r. In the main part of the paper one needs 
only the z component of the vector D^, i.e., Dzz- 
Using the integral formula 



dV rotF = dSxF 
V Js 



and the relation 



rot [/A] = [V/ X A] + /rotA, 
one can rewrite Eq. (jA6[) as 



[e^ X (r' - r)] 



(A8) 



(A9) 



D;(r) ^ f dS'x aziv') 
Js 



j3 /w_ „ X (r'-r)] 



dVVcr^(r') X 



r — r' I >ro 



r'-r 



rfS'xa.(r')^^^^^^^^^. (AlO) 



r' - r 



For a flat domain wall with 



dffziz') 



(All) 



the second term in Eq. (jA10[) is zero. In this case one 
has 

D;(r) = iy dS' X az{r') ^ \ ^ . - — a,(r)e„ 



JS |r - r| 

(A12) 

where the integral is taken over the surface of the sam- 
ple, the vector dS' is directed outwards, and the last 
term is the integral over the sphere ro, taken with 
(Jziv') => cTzir). The first term in Eq. (|A12p . multiplied 
by g^gS/vo, yields the macroscopic internal field in the 
sample that also can be obtained from the Biot-Savart 
formula with the molecular currents j =c rot M flowing 
on the sample's surface, M ={gngSv/vo)e.z inside the 
sample and M = outside the sample. 

In particular, for crz(r) — Uz — const for ellip- 
soidal samples the integral in Eq. (|A12p becomes 
47rj/ (l — n^^') cr^e^, where 71*^^^ is the demagnetizing 
factor. This yields D'^(r) = D^^^cr^e^, where D'^^ = 
47ri/ (1/3 - n("^)) . Then Eqs. and (|X5|) resuh in 

Eq. (|14p in the homogeneous case. For a cylinder mag- 
netized with Uz — (Jz{z), the field along the symmetry 
axis is given by Eq. ^TE\i . 
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